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Nanoscale inhomogeneity seems to be a central feature of the d-wave superconductivity in the 
cuprates. Such a feature can strongly affect the local density of states (LDOS) and the spectral 
weight functions. Within the Bogoliubov-de Gennes formalism we examine various inhomogeneous 
configurations of the superconducting order parameter to see which ones better agree with the ex- 
perimental data. Nanoscale large amplitude oscillations in the order parameter seem to fit the LDOS 
data for the underdoped cuprates. The one-particle spectral function for a general inhomogeneous 
configuration exhibits a coherent peak in the nodal direction. In contrast, the spectral function in 
the antinodal region is easily rendered incoherent by the inhomogeneity. This throws new light on 
the dichotomy between the nodal and antinodal quasiparticles in the underdoped cuprates. 

PACS numbers: 74.20.-z, 74.72.-h, 74.25.Bt 



I. INTRODUCTION 

The scanning tunneling microscopy (STM) and the 
angle-resolved photoemission spectroscopy (ARPES) are 
two of the most important tools for unraveling the mys- 
tery of the high-temperature superconductors (HTS). 
The existence of a Fermi surface and the d-wave sym- 
metry of the superconducting state are very important 
properties revealed by ARPES Q in momentum space. 
In contrast, STM has provided important complemen- 
tary information in real space through the measurement 
of the local density of states (LDOS). A surprising fea- 
ture of the HTS seen through STM is the conspicuous 
inhomogeneity EILU 

While some of the STM and ARPES data are straight- 
forward to interpret, others are not. The inhomo- 
geneity underlies much of the difficulty as there has 
not been much theoretical work addressing the effect 
of inhomogeneous d-wave superconductivity (DSC) on 
the LDOS IE 0, and the one-particle spectral func- 
tion HHEJ This 

paper is intended to partially remedy 

the situation. 

In the absence of a complete theory of HTS, what we 



have done is to examine various types of inhomogeneity 
for comparison with the STM and the ARPES data. In 
this way we hope to extract as much information from the 
data as possible. In the following, we first introduce the 
model Hamiltonian and the method of calculation. The 
calculated results for LDOS and the spectral function are 
presented and their implications are discussed. 



II. MODEL HAMILTONIAN 

In this work, we focus exclusively on the effect of inho- 
mogeneous d-wave pairing field. We therefore adopt the 
following Hamiltonian 

H = ^(e fe - ^)c{ a c ka + ( A *j c n c jl + H - c ) (1) 

ka <i,3> 

where cj CT creates an electron on site i with spin <r, and 
ct CT an electron with momentum k and spin a. < i,j > 
is a nearest neighbor pair. The kinetic energy is give by 



f-k = ti(cosk x +cosk y )/2 + t2Cosk x cosk y + t^,{cos2k x + cos2k y )/2 + t^{cos2k x cosk y + cosk x cos2k y )/2 + t<^cos2k x cos2k, 

I 



where the hopping parameters ti_5=-0.5951, 0.1636, 
-0.0519, -0.1117, 0.0510 eV are from Bi 2 Sr 2 CaCu 2 8+ s 
band structure determined by Norman et al. Ajj 
is the d-wave pairing amplitude over a nearest neighbor 
pair < i, j >. 

For any spatial distribution of A^- , the Hamiltonian (1) 



can in principle be straightforwardly diagonalized after 
performing a Bogoliubov transformation. In reality, due 
to the large lattice size involved we need to resort to a 
special technique for calculating the LDOS. 

Gagliano and Balseiro |l2( have proposed an efficient 
method for calculating the resolvent Ga —< "00 1 A'(Z — 
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H) 1 ^4|-0n >■ Ga(Z) is expressible as a continued frac- 
tion QQ 

G A (Z) = <^ A ' A ^> . (2) 
Z - a a J— p— 

where the coefficients a* and bi can be obtained from 
A \ipo > by repeated application of the Hamiltonian 
if. The Fourier transform of the self-correlation func- 
tion Ca{P - t) =< i)a\A^(t)A(t')\ipo > can be recov- 
ered from the imaginary part of the resolvent Ca (w) = 
iJmGifw + ijj + Bo). 

The above method is originally devised for calculat- 
ing the dynamical properties of quantum many-body sys- 
tems |l2l . but it can be easily adapted for our purpose. 
To calculate the LDOS, we simply take \ipo > to be the 
vacuum and A be c\ a . For the spectral function, we 
choose c\ a instead. The imaginary part of Ga then yields 
the LDOS and the spectral function respectively. 

III. LDOS 

For a homogeneous superconductor, the gap parame- 
ter can be directly determined from the measured LDOS. 
For an inhomogeneous superconductor, it is not triv- 
ial to invert the LDOS data for A.^. Our approach 
is to try various Ay configurations to fit the LDOS 
data. The first set of data we want to fit are the LDOS 
spectra measured at representative points in underdoped 
cuprates as shown in Figure 1. Motivated by the 
experimental gapmap, we consider a cone-shaped distri- 
bution Q of the A order parameter described by the inset 
of Figure 2. A rises to 4A at the center of a 400x400 
lattice, but it returns essentially to the background value 
Ao = 0.028eV for distances larger than ten lattice spac- 
ings. The calculated LDOS spectra at a series of points 
(at a distance 80,8,6,5,4,2,1 from the origin) are displayed 
in Figure 2. They indeed resemble the measured spectra 
in Figure 1. In particular, we see that the increase in 
gapsize as one moves toward the center of the lattice is 
accompanied by a gradual degradation of the coherence 
peak. As a reference, the LDOS of a uniform supercon- 
ductor with A = 4Ao is included in Figure 2. 

One implication of the above result is that the large- 
gap incoherent LDOS spectra characteristic of the un- 
derdoped cuprates can be interpreted as a rapid rise of 
the pairing field in a small region, or a nanoscale large 
amplitude fluctuation. There is no need to invoke a cer- 
tain unknown charge-ordered zero-temperature pseudo- 
gap state 0| which competes with DSC. 

The other notable feature of Figure 2 (which agrees 
with Figure 1) is that the low energy portions of the 
LDOS spectral are nearly identical suggesting homoge- 



neous nodal superconductivity coexisting with the inho- 
mogeneous antinodal feature. Such a contrast between 
nodal and antinodal excitations is also reflected in the 
ARPES spectra to be discussed in the next section. 

Besides the cone-shaped distribution in Figure 2, we 
have also considered many other A configurations which 
do not fit the data so well. Figure 3 shows the LDOS 
spectra corresponding to a mesa-like A configuration de- 
scribed in the inset (A = 4Ao within a distance of three 
lattice spacings away from the origin and A = Ao else- 
where). Such an extended region of high A leads to a 
higher energy coherent peak in the LDOS spectra near 
the center of the mesa. It is quite clear that the data in 
Figure 1 can discriminate such a configuration from the 
previous configuration of Figure 2. While the fit to the 
data in Figure 1 may not be unique, Figure 2 is the best 
one we have come up with so far. 

Although we have considered only a single cone-shaped 
A configuration, we anticipate that in real cuprates there 
would be a disordered array of cones. As long as the cones 
do not overlap strongly, the LDOS of the system should 
resemble that of a single isolated cone. 

The second set of data to fit are the Fourier- 
Transformed (FT) LDOS. Due to the constraint of com- 
puter time, we again consider only one cone. The cal- 
culated FT-LDOS spectra at various energies are dis- 
played in Figure 4 together with the measured one [lij . 
As emphasized by Dell' Anna et al. 0, previous analy- 
ses involving impurity scattering of quasiparticles tend 
to yield LDOS patterns with extended curve-like fea- 
tures 0, 0, 0, H3| in the high intensity regions in 
momentum space. This is in contrast to the spot-like 
intensity patterns seen experimentally. A zeroth order 
approximate calculation by DellAnna et al. shows that 
a mesoscopically inhomogeneous A distribution indeed 
yields a central spot. Our result in Figure 4 is exact and 
gives more details than theirs. In particular, the result 
exhibits high intensity features along the diagonal as well 
as horizontal and vertical directions in good agreement 
with experiment. The checkerboard-like intensity modu- 
lations of periodicity about four lattice spacings are also 
reproduced. 

IV. SPECTRAL FUNCTION 

The spatial inhomogeneity revealed by the STM data 
seems to be at odds with the well-defined Fermi arc. 
In addition, the nodal quasiparticle peak remains well- 
resolved even in strongly underdoped cuprates. In con- 
trast, the antinodal quasiparticle peak is well-defined 
only near optimal composition [ij, l2l| . This contrast has 
led to the speculation that the nodal and antinodal exci- 
tations have different origins 0] . One is associated with 
DSC which dominates in the optimal and overdoped re- 
gions, whereas the other is related to an unknown pseu- 
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dogap state which competes with DSC and which domi- 
nates in the strongly underdoped samples. Other expla- 
nations include coupling of the electrons with the (it, it) 
magnetic excitations |22( and a scattering mechanism op- 
erating mainly on the antinodal quasiparticles |2l) . 

We have seen in the previous section that the inco- 
herent LDOS spectra in the underdoped samples can be 
explained in terms of inhomogeneous DSC. Here we at- 
tempt to do the same for the ARPES spectra. Figure 5 
is the calculated spectral density for a disordered (ran- 
domly positioned) array of cones near the nodal direction 
(5a) and antinodal direction (5b), both momenta are lo- 
cated on the Fermi surface. The nodal quasiparticle peak 
is indeed well-resolved. The antinodal peaks are broader, 
but they are still resolved. This is because the A distri- 
bution in Figure 5 is dominated by low A values. 

To simulate the incoherent antinodal spectral func- 
tions seen experimentally in underdoped samples, we fab- 
ricate a more disordered A configuration in Figure 6. 
The antinodal spectra are indeed incoherent, whereas the 
nodal one remains coherent. We have examined other A 
configurations, the contrast between the nodal and antin- 
odal spectral features seems to be generic for inhomoge- 
neous DSC independent of the details of the inhomo- 
geneity. Such a result is actually reasonable because the 
nodal quasiparticles have vanishingly small excitation en- 
ergies independent of the magnitude of A, therefore they 
can propagate freely in any inhomogeneous superconduc- 
tor. For the opposite reason, the energy of an antinodal 
quasiparticle is very sensitive to A, high energy antinodal 
quasiparticles are confined to regions of high A. 

V. DISCUSSION AND CONCLUSION 

For simplicity , we have limited ourselves to only one 
kind of inhomogeneity in this paper, the inhomogeneous 
pairing field. In real cuprates, impurities are present as 
well as short-range antiferromagnetism. They could also 
affect the LDOS 0, an d the spectral weight function. 
Further theoretical study is required to include those ef- 
fects. 

Despite the limitation, our results so far support the 
following conclusions: (1) Inhomogeneous DSC is an im- 
portant determining factor in the LDOS and ARPES 
spectra; (2) It can explain many unusual features of the 
experimental spectra without invoking an unknown state 
which competes with DSC, at least for optimally doped 
and moderately underdoped systems. The "dichotomy" 
between nodal and antinodal excitations seems to be a 
mere consequence of the inhomogeneous d-wave pairing 
field. 
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Fig 1: LDOS spectra measured at representative points 
in underdoped cuprates, data taken from ref. [lij . 

Fig 2: Calculated LDOS spectra for the spiky A con- 
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figuration described in the inset. The curve labeled 'uni- 
form A' is included as a reference. It is the LDOS spec- 
trum for a uniform A = 4A . 

Fig 3: Calculated LDOS spectra(at a distance 
43,4,3,2,1,0 from the origin) for the A configuration de- 
scribed in the inset. 

Fig 4: Calculated FT-LDOS (left column) for a 



cone-shaped A distribution at various biases (w — 
—8, —16, —26 meV from top to bottom) compared with 
data (right column) taken from ref. |15j . 

Fig 5: Spectral functions along the nodal and antinodal 
directions corresponding to the A configuration on the 
left panel. 

Fig 6: Same as Figure 5 for a different A distribution. 
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This figure "fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0508358vl 



o 


30 - 


o 


25 - 






=3 
LL 


20 - 




15 - 




i_ 


10 - 


O 


0) 


5 - 


Q. 


<J) 


- 



(a) 



(0.397t , 0.3971) 



(b) 



(it , 0.257t) 



-3-2-1 1 2 3-4 -3 -2-10 1 

Binding Energy (100 meV) 



2 3 4 



This figure "fig6a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0508358vl 
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